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ABSTRACT 


The  purpose  of  this  discussion  is  to  relate  the  previously  developed 
matrix  analysis  of  electromagnetic  scattering  explicitly  to  the  experi¬ 
mentally  measurable  bistatic  polarization  matrix.  In  addition,  in  order  to 
perform  the  numerical  computations  more  conveniently,  a  new  matrix  analy¬ 
sis  is  described,  in  which  the  constraints  of  reciprocity  and  energy  conser¬ 
vation  are  employed  in  the  course  of  obtaining  a  solution. 


iii 


ACKNOWLEDGEMENT 


The  author  is  indebted  to  several  colleagues:  J.  F.  A.  Ormsby, 
W.  T.  Payne,  M.  R.  Weiss,  L.  J.  Kaplan  and  C.  V.  McCarthy  for 
many  helpful  discussions  during  the  course  of  this  work. 


iv 


TABLE  OF  CONTENTS 

SECTION  I 

Page 

INTRODUCTION  1 

SECTION  II 

THE  BISTATIC  POLARIZATION  MATRIX  3 

SECTION  III 

SOLUTION  OF  THE  SYSTEM  OF  MATRIX 

EQUATIONS  12 

REFERENCES 

17 

V 


SECTION  I 


INTRODUCTION 


The  purpose  of  this  discussion  is  as  follows: 

(1)  To  relate  the  previously  developed  matrix  analysis  of 
electromagnetic  scattering  explicitly  to  the  experi¬ 
mentally  measurable  bistatic  polarization  matrix. 

(2)  To  describe  a  new  matrix  analysis,  suitable  for  numerical 
computation,  in  which  the  constraints  of  reciprocity 

and  energy  conservation  are  employed  in  the  course  of 
obtaining  a  solution. 

Description  of  the  measurable  far  field  quantities  requires 
mapping  a  vector,  the  polarization  of  the  incident  plane  wave,  into 
another  vector,  giving  the  polarization  of  the  scattered  field; 
this  operation  is  conveniently  described  by  the  dyadic  amplitude  £  . 
In  Section  II  the  spatial  electromagnetic  fields  are  first  set  up  in 
terms  of  a  theoretical  quantity,  the  transition  matrix  T  ,  then 
the  dyadic  amplitude  is  defined,  and  by  specifying  bases,  e.g., 
linear-linear,  the  experimentally  measurable  bistatic  polarization 
matrix  P  is  obtained.  Finally,  the  determination  of  the  dyadic 
amplitude  from  the  transition  matrix  is  described. 

Section  III  gives  the  matrix  analysis,  in  which  contraints  of 
reciprocity  and  energy  conservation  are  incorporated  into  the 
theoretical/numerical  determination  of  the  transition  matrix  T  . 
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In  practice,  this  approach  may  make  feasible  the  employment  of  the 
existing  MITRE  computer  program  for  a  wider  range  of  target  sizes  and 
shapes  than  heretofore  possible. 

It  should  perhaps  be  pointed  out  that  the  discussion  throughout 

is  expected  to  be  applicable  to  general  radar  targets  with  no  energy 

dissipation  in  otherwise  free  space,  although  to  date  the  requisite  Q 

matrix  of  Section  III  has  only  been  discussed  for  general  conducting 

bodies  and  has  only  been  implemented  on  the  computer  for  conducting 
bodies  having  an  axis  of  rotational  symmetry. 
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SECTION  II 


THE  BISTATIC  POLARIZATION  MATRIX 


We  work  in  the  matrix  formalism  described  earlier  L^Jusing  the 
recent  modifications.  A  more  detailed  description  of  some  of  the 

quantities  discussed  below  may  be  found  in  the  references. 


An  arbitrary  monochromatic  incident  electromagnetic  field  (time 
dependence  exp  (-icot)  suppressed)  with  no  singularities  in  the 
immediate  vicinity  of  a  conducting  target  may  be  described  by  the 
electric  field  vector  (underline  denotes  vector  quantity) 


E1(r) 


V 

L 

omn 


£  (2n  +  1)  (n  -  m)  ! 
m 


n(n  +  1) (n  +  m) ! 


-1/2 


a  -  (r)  +  b  -  (r) 

omn  omn  —  omn  omn  — 


(1) 


at  any  field  point  _r  measured  from  an  origin  lying  in  the  interior 

of  the  target.  The  a  and  b  are  complex  expansion  coefficients, 

°  omn  omn  1 


£  is  the  Neumann  factor  e  =  1 ,  £  =2  otherwise,  and  the 

m  o  m 


spherical  partial  waves  are  given  explicitly  by 


[2] 


M  =  M  =  Vx 
”umn  — emn 
o 


r^  cos  mcf)  Pm  (cos  0)  j  (kr) 
sin 


=  [n(n  +  1) ] 


1/2 


£e»n  <9'W  Vkr) 

o 
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N_  =  N  =  (1/k)  V  x  M 

~ Cmn  ~emn  — emn 

o  o 


n(n  +  1)  P.  emn  (0,4>)  jn  (kr)  /  (kr) 
o 


(2) 


+  [n(n  +  l)]1^2  (0,40  (1/kr)  d[rjn(kr)  ]/dr 

o 


The  three  indices  run  through  values  o  =  e,  o  (even  or  odd), 

m=0,  1,  ...  ,  n,  and  n  =  1,  2,  ...  .  The  Pm  are  associated 

n 

Legendre  functions,  and  the  jn(kr)  are  spherical  Bessel  functions. 

r  3  1 

The  vector  spherical  harmonics  are  defined  byL  [note  that 
r  =  r(0  ,(j>)  ] 


P  (r) 
—emn 
o 


%  cos  -r^m  / 
r  m<i>  P  (cos  0) 

sin  n 


B  (r) 
-emn 
o 


[n(n  +  1)]  r  V[C°S  racj>Pra  (cos  0)  ] 

sin  n 


f  x  C 

-emn 

o 


(3) 


C  (r) 
—emn 
o 


[n(n  +  1)  ]  1//2  V  x  [_r 

—  sin 


Pm  (cos  0)] 
n 


-  r  x  B 

—emn 
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In  similar  fashion,  the  most  general  scattered  wave  satisfying  the 
radiation  condition  at  infinity  may  be  written 


IS(r)  =  l  i 
omn 


e  (2n  +  1)  (n  -  m)  ! 
m _ 

n(n  +  1) (n  +  m) ! 


1/2 


fomn  +  8-™ 


J<3) 

’omn 
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(4) 


in  terms  of  the  outgoing  spherical  vector  partial  waves,  obtained  by 

replacing  Jn(kr)  by  h^kr)  (the  spherical  Hankel  function  of  the 

first  kind)  everywhere  in  Equation  (2). 

Because  of  the  linearity  of  the  problem,  the  coefficients 

f  ,  g  of  the  scattered  wave  must  be  obtainable  by  a  linear 

omn  omn 

operation  on  those  of  the  incident  wave:  one  writes 


f  -  l  [T1  ,  ,  ,  a  ,  ,  ,  +  T2  ,  ,  ,  b  ,  ,  ,  ] 

omn  .  v  .  omno  m  n  a  m  n  amna  m  n  a  m  n 

a  m  n 


(5a) 
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g  -  "  [Td  f  f  ,  a  f  f  f  +  T4  ,  ff  t  b  t  «  t 1 

Gomn  omno  m  n'  o  m  n1  omno  m  n'  o  m  n1 

a’m'n* 


or,  in  an  obvious  matrix  notation 


f 

g 

The  matrices  having 
where  the  indices  a,  m, 
or  more  briefly 

T 


pi  p2 

a 

p3  p4 

b 

elements  ^  \2 

n  have  been  reordered 

p  1  p2 

p3  p4 


(5b) 

,  and  so  forth, 
into  a  single  index, 

(5c) 
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will  be  called  VlCLnbxdXOn  m&Oiix  .  Once  obtained  for  a  given  target, 
it  enables  one  to  compute  the  scattered  wave  corresponding  to  any 
given  incident  wave  (1) ,  by  applying  Equations  (4)  and  (5) . 

In  radar  work  one  is  primarily  interested  in  the  scattered  far 
field  due  to  excitation  by  a  plane  wave  having  specified  direction 
of  incidence  and  specified  polarization.  Because  we  are  interested 
in  mapping  a  vector  (the  incident  polarization)  into  another  vector 
(the  far  scattered  E  field)  it  is  convenient  to  introduce  the 

[4  51  - 

dyadic.  L  5  J  far  field  amplitude  k^)  describing  the  scattered 

electric  vector  far  from  the  target  (i.e.,  kr  >>  1).  The  directions 
of  incidence  and  scattering  (i.e.,  observation)  are  given  by 

A  A  A  ^ 

k.  55  k.(0.,  $  )  and  k  =  k(0  ,  4>)  »  respectively,  and  e.  denotes 
1  XX  1  S  S  S  S  1 

the  electric  polarization  vector  associated  with  the  incident  wave. 

One  then  can  write 

E^(_r)  ^  (1/kr)  exp(fkr)  F  (kg,k^)  ’  as  kr  00  ,  (6) 

Because  both  the  incident  wave  and  the  far  field  E  are  transverse 
to  their  directions  of  propagation,  the  dyadic  amplitude  may  be 
constrained  to  satisfy 

F(k  ,k. )  *  k  -  k  •  F(k  ,k.)  =  0  (7) 

=  s  i  i  s  =  s  i 

thus  having  only  four  independent  components •  In  addition,  the 
reciprocity  theorem,  relating  to  interchange  of  source  and  observation 
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point,  requires  that 


A  A  T*  A  A 

f(ks,k.)  =  F  (-k.,  -kg) 


(8) 


where  the  superscript  T  designates  the  dyadic  transpose.  It  may 
also  be  shown,  as  a  consequence  of  energy  conservation,  that^ 


dfi  F1  (k  k  )  '  F(k  ,k  )  =  -2tti  [F(k  ,k'  )  -  fT  (k *  k  )] 
s—  si  —  si  —  11  —  11 


(9) 


where  the  integration  is  over  the  unit  sphere  of  directions  of  . 


The  relationship  between  the  dyadic  amplitude  and  its  constraints, 
on  the  one  hand,  and  the  transition  matrix,  on  the  other  hand,  will 
be  made  apparent  shortly.  It  is  convenient  first,  however,  to 
examine  the  (2x2  complex)  bistatic  polarization  matrix.  One  can 
define  the  vector  scattering  amplitude  e^  to  be  just  ( r) 

exclusive  of  the  radial  factor  (1/kr)  exp(ikr)  ,  so  that  from  Equa¬ 
tion  (6)  there  results 


e  =  F(k  ,k.)  *  e. 
— s  =  s  i  — 1 


(10) 


The  incident  polarization  may  be  written 


e . 


i 


(ID 


i 


7 


where  if  both  components  are  real  (or  if  both  are  complex,  but  with 
equal  phase)  the  incident  wave  is  linearly  polarized.  Otherwise,  the 
incident  wave  is  elliptically  polarized. 

In  similar  fashion  the  scattered  polarization  vector  may  be 

written 


e  =  e_  0  +  e  <j> 

“S  0S  d>  Ts 
s  Ys 


(12) 


and  this  of  course  will  in  general  represent  an  elliptic  state  of 
polarization,  even  though  the  incident  wave  may  be  linearly  polarized. 
The  base  vectors  of  Equations  (11)  and  (12)  appear  the  natural  choice 
for  a  linear-linear  basis,  in  terms  of  which  the  2x2  bistatic 
polarization  matrix  P  enables  one  to  compute 


e^ 

0 

s 

e 

d> 

Ys 

P  P 

0  9  0  cp 

si  s 


\  0t  \ 


'0. 

l 


(13a) 


in  which 


P  q  =  0  •  F  •  0. 

0  0.  s  i 

s  i 


P  A  _  0  •  F  •  A. 

0  cp .  =  S  =  V1 

s  1 


(13b) 


P  =  (J)  •  F  *  6 . 

<^s  9i  8=1 


P  =  d>  •  F  •  <b. 

4  cp.  vs 

Ys  1 
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Two  observations  should  be  noted  in  the  system  of  Equations 
(10  -  13).  First,  the  dyadic  amplitude  F  is  in  a  sense  more  general 
than,  say,  the  P  of  Equation  (13),  in  that  F  provides  a  full  des¬ 
cription  of  the  problem  without  having  specified  the  bases,  e.g., 
linear-linear,  circular-circular,  linear-circular.  The  other  obser¬ 
vation  is  that,  in  considering  the  bistatic  problem,  one  is  led  quite 
naturally  to  choose  two  distinct  coordinate  systems  for  the  description 
of  incident  and  scattered  polarization  states.  This  has  the  notable 
consequence  that  when  one  specializes  P  to  the  monostatic  case, 
choosing  k  =  -  k.  (for  which  one  easily  verifies  that  0  =  0.  , 

4>g  =  -<fK  )  the  linear-linear  polarization  matrix,  say  for  a  conducting 
sphere,  becomes  (  |a|2  =  radar  cross  section) 


P  = 


a 

0 


0 


-a 


(14) 


This  result  differs  from  the  accepted  (single  coordinate  system)  des¬ 
cription  of  monostatic  scattering,  in  which  the  minus  sign  is  absent. 
Regardless  of  accepted  convention,  however,  the  two-coordinate  system 
description  giving  rise  to  Equation  (14)  appears  almost  imperative 
wherever  bistatic  as  well  as  monostatic  measurements  are  being 
considered . 

So  far  we  have  seen  how  the  dyadic  amplitude  F  (k^  ,  k^), 
which  might  be  regarded  as  the  basic  quantity,  relates  to  the  bis¬ 
tatic  polarization  matrix  P  and  hence  experimental  measurement.  It 
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only  remains  to  relate  F  to  the  theoretical  formulation  of  the 
problem,  in  particular  to  the  transition  matrix  T  .  The  desired 
result  is  obtained  almost  immediately,  by  1)  introducing  the 
appropriate  plane  wave  expansion  coefficients  ^  in  the  incident  wave 
of  Equation  (1,2)  employing  the  asymptotic  form  of  the  spherical 
Hankel  functions  valid  at  large  distances  from  the  origin  ^  ^  in  the 
scattered  wave  of  Equations  (4),  and  finally,  (3)  comparing  the  result 
with  Equation  (6)  to  get 


F(ks  ,  k  )  =  l  <-i) 


.  \  n 


r 


Omn 
a'm'n' 


e  (2n  +  1)  (n  -  m)  ! 
m 


(n  +  m)  ! 


1/2 


(15) 


•p  1  p2 

C(k±) 

iB(ks)] 

p3 

-iB(k.) 

d)n' 


e  ,(2n '  +  1)  (n'  -  m’)  ! 
m  1 _ _ _ _ 

(n ’  +  m! ) ! 


where  summation  indices  have  been  suppressed  in  the  expression  in  curly 
brackets  for  clarity.  This  expression  when  written  out  consists  of 
the  sum  of  dyads  of  four  types,  the  first  of  which  is  given  by 


C 

“omn 


(k  )  T 
s 


1 

cmna ’ mf n 1 


“”tj  ’mfn 


and  so  forth. 

Because  the  spherical  vector  harmonics  C(k  )  and  B(k  )  have 

—  s  —  s 

no  components  in  the  observation  direction  ,  the  constraints  of 
Equation  (7)  are  seen  to  be  satisfied  by  inspection.  Comparison  of 


'  1/2 


10 


Equation  (15)  with  the  reciprocity  theorem  Equation  (8) ,  however,  in 

[3j 

view  of  orthogonality  of  the  spherical  vector  harmonics  1  J  ,  leads  to 
the  requirement  that  the  transition  matrix  be  symmetric,  i.e.. 


T  *  T* 


(16) 


or,  in  more  detail, 
Finally,  the  energy 
constraints  that  T 


(T1)'  =  T1  ,  (T2)'  =  T3  ,  (T3) '  =  T2  ,  (T4) '  = 
conservation  requirement  (9)  gives  rise  to  the 
be  unitary-related,  i.e.,  that 


T4. 


T'*T  =  Re  T  . 


(17) 
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SECTION  III 


SOLUTION  OF  THE  SYSTEM  OF  MATRIX  EQUATIONS 

Rather  than  deal  with  the  transition  matrix  T  ,  the  procedure 
to  be  described  is  slightly  more  transparent  if  one  employs  the  S- 
matrix  defined  by 


S  =  1  -  2T  . 


(18) 


The  derivation  of  a  matrix  equation  for  the  determination  of  S  has 
been  discussed  elsewhere.  Specifically,  one  has 


QS  -  -Q* 


(19) 


where  the  matrix  elements  of  Q  are  explicitly  known.  To  this 
equation  may  be  adjoined  the  constraints  due  to  Equations  (16)  and 
(17) ,  which  are  readily  seen  to  be 


S  *  Sf  (20) 

and 

S'*S  -  1  (21) 

i.e.,  S  is  symmetric  and  unitary,  respectively*  Two  extremes  of  view 
with  regard  to  the  system  of  Equations  (19)  ,  (20)  and  (21)  are  as 
follows:  First,  one  might  truncate  the  infinite  matrix  Equation  (19), 
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solve  numerically  by  digital  computer,  then  compare  the  resulting 
solution  with  Equations  (20)  and  (21) ,  the  latter  thus  being  employed 
as  consistency  checks.  On  the  other  hand,  one  might  attempt  to  treat 
all  three  equations  from  a  unified  point  of  view  from  the  onset, 
obtaining  a  solution  in  some  sense  of  Equation  (19)  subject  to  the 
constraints  (20)  and  (21) .  The  first  approach  has  been  employed  up  to 
now  on  the  computer  for  bodies  of  rotational  symmetry,  and  works  quite 
satisfactorily  for  a  restricted  range  of  body  shapes  and  sizes.  The 
second  approach  is  proposed  in  an  effort  to  extend  the  range  of 
bodies  that  can  be  handled,  in  view  of  the  fact  that  the  constraints 

2 

essentially  determine  three  quarters  of  the  solution  [i.e.,  of  the  2N 
real  parameters  appearing  in  the  N  x  N  (truncated)  complex  matrix  S  , 
it  can  be  shown  that  only  N(N  +  l)/2  are  independent,  if  S  satisfies 
( 20  and  21)  ] . 

To  develop  a  unified  analysis,  observe  first  that  if  S  could 
be  constructed  in  the  form 


S  =  U1  U 


(22) 


where  U  is  unitary,  then  both  constraints  would  be  satisfied  by 
inspection.  This  suggests  that,  rather  than  inverting  Q  directly 
in  Equation  (19),  it  be  made  unitary.  Thus,  consider  the  upper  tri¬ 
angular  matrix  M  (i.e.,  all  elements  are  zero  below  the  main  diagonal) 
which,  by  premultiplication,  makes  Q  into  a  unitary  matrix  Qun^t  > 
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vie  . 


MQ 


Q 


unit 


(23) 


Premultiplying  Equation  (19)  by  M  ,  one  can  write 


Q  .  S  =  -  MQ*  =  -MM*  1  Q*  . 
unit  x  x  unit 


Upon  solving  for  S  there  now  results 


S 


-Q'*  (MM*  1)  Q*  . 
unit  unit 


(24) 


(25) 


Substituting  this  result  in  Equation  (20),  the  symmetry  constraint,  it 
follows  without  difficulty  that  the  matrix  product  MM*  ^  must  be 
symmetric.  But  each  of  the  matrices  appearing  in  the  product  is 
upper  triangular,  and  their  product  is  again  upper  triangular.  Con¬ 
sequently,  the  product  must  be  a  diagonal  matrix.  Further,  the 
diagonal  elements  can  be  written  out  explicitly,  giving 


MM* 


M  /M  *  0 

ir  11 


0  M  /M  * 
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If  next  we  can  arrange  to  choose  the  diagonal  elements  of  M  to  be 
real,  then 


* 


MM* 


-1 


=  1 


(25) 


From  Equation  (24)  the  S-matrix  is  now  given  by 


S 


Q  1  *  Q* 

unit  unit 


(26) 


which  is  of  the  required  form  (22).  Substituting  Equation  (26), 

along  with  the  identity  Qf*  .  Q  .  =  1  back  in  Equation  (18), 

J  unit  unit 

the  desired  transition  matrix  is  finally  given  by 


T  -  Qf’ 


unit 


Re 


(<W 


(27) 


Returning  to  M  for  a  moment,  Equation  (25)  states  simply  that 
M  is  real.  Thus,  the  process  may  be  summed  up  in  the  (formal) 
theorem:  Given  the  matrix  Equation  (19)  ,  with  constraints  (20)  and 

(21)  on  the  solution,  it  follows  that  the  given  matrix  Q  cannot 
be  arbitrary,  but  must  be  such  as  to  be  factorizable  into  the  product 
of  a  real  upper  triangular  matrix  and  a  unitary  matrix,  namely 


Q  = 


Q 


unit 


(28) 
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In  order  to  carry  out  numerically  the  above  analysis,  which  is 
formal  in  the  sense  that  infinite  matrices  have  been  treated  without 
rigorous  verification  of  the  consequent  limiting  processes  involved, 
it  is  of  course  necessary  to  truncate  S  and  Q  at  a  finite  number 
N  of  rows  and  columns.  Having  done  this,  construction  of  the 
N  x  N  unitary  matrix  Qun^t  is  done  by  Schmidt  orthogonalizat ion 
of  the  N  vectors  given  by  the  rows  of  Q  ,  beginning  with  the  bottom 
row  and  working  up.  Thus,  writing  for  the  N  component  vector 

obtained  from  the  bottom  row,  is  first  normalized  to  unit  length 

by  multiplying  by 


(V  •  V 


-1/2 


^*Nn  ^Nn 


-1/2 


This  constant,  which  is  real  by  definition,  is  just  the  last  diagonal 

entry  of  the  triangular  matrix,  whereas  the  normalized  vector 

constitutes  the  bottom  row  of  Q  Next  one  chooses  a  linear 

unit 

combination  of  and  the  vector  _  obtained  from  the  next 

to  last  row,  so  that  the  resulting  vector  is  orthogonal  to  and 

so  forth. 

A  brief  numerical  investigation  of  this  procedure  has  been  per¬ 
formed  in  order  to  ascertain  its  feasibility  and  usefulness  in  the 
computer  program  for  electromagnetic  scattering.  Results  indicate  more 
rapid  convergence,  versus  truncation  size  N,  than  is  obtained  by 
straightforward  numerical  solution  of  Equation  (19)  employing  matrix 
inversion . 


0 
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